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ABSTRACT 

We have calculated the 2.2 MeV gamma-ray line intensity from the 

Sun using a Monte-Carlo method for neutron propagation in the solar 

atmosphere. We provide detailed results on the total gamma-ray yield 

per neutron and on the time profile of the 2.2 MeV line from an 

instantaneous and monoenergetic neutron source. The parameters which 

have the most significant effects on the line intensity are the energies 

of the neutrons, the position of the neutron source on the Sun, and the 
3 

abundance of He in the photosphere. For an xsotropic neutron source 

which is not too close to the limb of the Sun, the gamma-ray yield is 

between about 0.02 to 0.2 photons per neutron, provided that the 

3 

neutron energies are in the range 1 to 100 MeV and the ratio He /H is 

_5 

less than about 5x10 
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I- Introduction 

Gamma-ray line emissions at 0.51, 2.2, 4.4, and 6.1 MeV were 

observed by Chupp et_ al^ (1973) during the flare of August 4, 1972. 

These emissions are believed to be due to positron annihilation, 

deuterium formation by neutron capture on hydrogen, and the deexcita- 

12 16 

tion of excited states in C and 0 (Lingenfelter and Ramaty 
1967, Ramaty and Lingenfelter 1973), 

The calculation of the line intensities from excited nuclei is 
in general straightforward because the deexcitations occur in a much 
shorter time scale than any of the characteristic times of the solar 
flare. However, because the time scales associated with positron 
annihilation and neutron capture in the solar atmosphere could be 
comparable with the characteristic times of the flare, the calculation 
of line intensities at 0.51 and 2.2 MeV has to take into account in 
detail the production, slowing down, and annihilation or capture of 
positrons or neutrons. In this paper we present a calculation on the 
formation of the 2.2 MeV line from neutrons in the solar atmosphere. 

The principal modes of neutron production in solar flares are 
nuclear reactions of accelerated particles with the ambient solar 
atmosphere (Lingenfelter ejt al. 1965) . Typically, the neutrons have 
energies between 1 and 100 MeV, but their number and energy spectrum 
depend on the intensity, composition, and the spectrum of the accelerated 
charged particles and the composition of the ambient medium in the flare 
region. Some of the effects of these variables were investigated by 
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Lingenfelter and Ramaty (1967) and Ramaty and Lingenfelter (1973) . 

However, in calculating the gamma-ray yield at 2.2 MeV, these authors have 
assumed that all downward moving neutrons are captured by protons, 
that all the upward moving photons escape from the Sun, and they have 

3 

neglected neutron capture on photospheric He . 

In this paper we present corrections to these simplifying assump- 
tions . We find that some of the downward moving neutrons decay or 

3 

escape before capture, and that He can constitute an important non- 
radiative sink for the neutrons. Furthermore, some of the upward 
moving photons are lost from the 2.2 MeV line by Compton scattering. 

We also calculate the time profile of the 2.2 MeV line intensity from 
an impulsive release of neutrons at the Sun, and we show how these 
calculations can be used to find the time -dependent 2.2 MeV line 
intensity from an arbitrary solar neutron source. 

II . Method of Calculation 

We consider the instantaneous release of monoenergetic neutrons 
with given angular distribution at a given height, h, above the photo- 
sphere. By using a Monte-Carlo method, we follow the path of each 
neutron in the solar atmosphere. This path is terminated by decay, 
escape or capture. For those neutrons which get captured by protons, 
we record the elapsed time from release to capture, and the depth in 
the solar atmosphere where the capture takes place. Thus, we can 
calculate the probabilities for decay, escape, and capture, the 
attenuation of the gamma rays, and the time profile of the gamma-ray 
burst. The details of the calculation are given in the Appendix. Here 
we only present our assumptions and approximations. The results are 
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presented in the next section. 

We assume that the neutrons are produced in the chromosphere or 
lower corona where the flare accelerated particles are expected to 
interact with the solar atmosphere. Then the value of h is small 
compared to the radius of the Sun, and hence the solar atmosphere can 
be treated as a plane stratified medium. However, because the exact 
value of h is not known, we treat it as a free parameter. 

Both the temperature and density in the atmosphere are functions 
of height, and we have used the tabulation of Allen (1963) for these 
quantities. The principal interactions between the neutrons and the 
solar atmosphere are elastic scattering and capture. Because of the 
relative abundances of the constituents of the solar atmosphere and 
their cross sections, elastic scattering is mainly due to ambient 
protons, but neutrons are captured on both protons and helium-3. The 
cross section, cr s , for elastic scattering is essentially constant at 
2 x 10 cm from 1 eV to 10 eV 5 and then it drops to about 10~ 23 cm 2 
at 100 MeV (Hughes and Schwartz 1958) „ The cross section for capture 
on protons is inversely proportional to velocity and is given by 
cr c (H)^2.2xl0 3 °6 1 cm 2 ? where cp is the relative velocity between the 
neutron and the proton. The capture cross section on He 3 has the same 
form as that on H, and is given by cr c (He 3 ) = 3.7 x 10 -26 B' 1 cm 2 (King 
and Goldstein 1949). 

Only upper limits are available on the He 3 abundance in the photo- 
sphere. In the solar wind, the ratio He 3 /H is about 5 x 10‘ 5 (Geiss 
and Reeves 1972). In our calculation we use as limiting values 
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He 3 /H - 5 x 10‘ 5 and He 3 /H = 0. 

In neutron-proton elastic scattering we assume a uniform angular 
distribution for the recoil neutron in the center of mass system* This 
is valid for neutrons of energies less than 10 MeV* In the 10 to 100 
MeV range, the distribution is approximately symmetric about 90°, and 
the assumption of uniform distribution does not effect significantly 
the final probabilities* 

In calculating the gamma-ray line intensity at Earth from the 

reaction n -f-p^d+7 at the Sun, we assume that the gamma rays are produced 

isotropically* The principal interaction between 2.2 MeV photons and 

the solar atmosphere is Compton scattering* From the Klein-Nishina 

-25 2 

formula, cr (2.2 MeV) = 1.4 x 1*0 cm * 

KN 

III. Results and Discussion 

The probabilities for decay, escape, and capture on protons and 

3 3 

loss on He are shown as solid lines in Figures 1 and 2 for He /H = 

-5 3 

5 x 10 and He /H = 0, respectively . In these calculations the neutrons 

9 

were released at h - 3 x 10 cm above the photosphere with a uniform 
angular distribution. Because the calculations are insensitive to the 
exact value of h as long as it is much smaller than one solar radius, 
we do not show results for other values of h. 

The capture and loss probabilities increase with increasing energy, 
because at higher energies the neutrons penetrate deeper into the 
photosphere. This prevents them from escaping and reduces the escape 
probability; it also leads to a shorter capture time, thereby reducing 
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3 -5 

the decay probability. When He /H = 5 x 10 , the probability for 

3 

loss on He almost equals the capture probability on protons. The 
escape probability is greater than 0.5, because all initially upward 
moving neutrons escape from the Sun. Note that the sum of all probabilities 
equals 1. 

Now we consider the gamma-ray intensities. Let sHt-t ,0 E ) 

o’ ’ n' 

be the time profile at Earth of the 2.2 MeV gamma-ray line from one 

neutron of energy E^ released at time t at the Sun measured in photons 
-2 -1 

cm sec , where 0 is the angle between the earth-sun line and the 

vertical to the plane stratified medium. We define the photon yield 

per neutron, f(0,E ), as 
n 

f ( e » E n ) =4 itR 2 J^ <i(t-t oS e s E n )dt , (1) 

O 

and we show it by the dashed lines in Figures 1 and 2. 

At low neutron energies and 0 near zero, f is close to the capture 

probability on protons. This means that gamma rays from low-energy 

neutrons observed close to the vertical escape essentially unattenuated 

from the Sun. At higher energies and at larger angles, however, there 

is significant attenuation of the gamma rays. 

In previous calculations, it was assumed that all downward moving 

neutrons are captured and that all upward moving photons escape from 

the Sun. In this case f should equal 0.5. However, from Figures 1 

and 2 we see that, depending on the energy of the neutrons, the location 

of the flare on the Sun (which determines the angle 9), and the amount 
3 

o He in the photosphere, this assumption can overestimate the gamma- 
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ray yield by at least a factor of 2,5, Furthermore, if the flare 
occurs close to the limb of the Sun, the 2.2 MeV line becomes essentially 
unobservable. For such limb flares, however, the 0.51 MeV line and the 
nuclear deexcitation lines will still be observable because these lines 

are believed to be produced above the photosphere. 

2 3 

The function 4ttR £ is shown in Figures 3 and 4 for He /H - 

-5 3 

5 x 10 and He /H = 0, respectively. We note that the rise times of 

the gamma-ray burst are less than about 25 seconds in all cases. These 

rise times increase with decreasing neutron energy and increase with 

increasing 9. We have prepared extensive tables of £ for a variety of 

3 

values for the parameters E , 0, and the ratio He /H. These tables 

are available in unpublished form and will be presented in a forthcoming 

dissertation. 

The function £ can now be used to calculate the time-dependent 

2.2 MeV line intensity at Earth, F(t,0), from a continuous non- 

monoenergetic neutron source at the Sun. Let q (t , E ) be the 

o’ n 7 

instantaneous rate of neutron production in the energy interval dE 

n 

around E^. Then F(t,Q) is given by 

F(t,6) = f° dE f* q(t ,E ) «4(t-t ,9 ,E )dt . (2) 

In the present paper we did not carry out a detailed study of F because 
the observational data on the time profile of the 2.2 MeV line is rather 
meager. Reppin et al. (1973), however, have estimated the time -dependent 
2.2 MeV intensity for the flare of August 4, 1972 by replacing fi(t-t ,9,E ) 
in equation 2 with the function (4ttR ) exp [-\(t-t )]. Here 
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X = T c + T d > where is the capture time of neutrons on hydrogen 

and Tj is the mean life time of the neutron. This approach neglects 

the dependence of <f> on 8 and E n> as well as the capture of neutrons on 

3 

He ; also, it does not take into account the Compton scattering of the 
2.2 MeV photons. Nonetheless, an exponential function can be used as 


S- crude approximation of provided that the capture time t is 

c 


approximately given by 


T c (sec) 


1.5 x 10 


where n is the density of the ambient solar atmosphere where most of 
the captures take place (Ramaty and Lirigenfelter 1973). These authors 
have argued that the bulk of the captures take place at a depth in the 
photosphere where the neutron suffers its first elastic scattering. 

Our calculations essentially substantiate this assumption and indicate 
that the captures take place at depths of a few hundreds kilometers 
below the photosphere . 
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APPENDIX 

We now describe the method for the evaluation of the probabilities 

for escape, decay or capture of a neutron that is produced in the 

solar atmosphere, and the calculation of the depth in the photosphere 

where the neutron is captured by a proton. We then use this depth to 

determine the probability for Compton scattering of an upward moving 

photon resulting from this capture. 

Consider a neutron of energy E 0 produced at a height h Q above the 

photosphere at time t e moving at a directional cosine 4 0 between its 

velocity vector and the vertical to the plane stratified solar atmosphere. 

(ju > o is defined as motion upwards) Let be the height at which 

the j th elastic scattering between the neutron and an atmospheric 

proton has occurred at time t so that after the scattering the neutron 

has energy E. and is moving with directional cosine u . • We want to find 
J J 

a new height h^ + ^ w ^ ere t * ie next elastic scattering could occur at time 

t^_^. We also want to determine whether the neutron decays or is 

3 

captured (on H or He ) between h^ and and t0 find the neutron 

energy and direction after scattering at h ^ provided it has not 
decayed or has not been captured. 

In the following we shall have to determine random numbers X with 
distribution functions $ (X) from uniformly distributed random numbers 
R, 0^R<1, which can be generated on a computer. In the subsequent 
discussion all R*s represent such random numbers. From the theory of 
probability 


$(X) = R and X = § - 1 (R) 


(Al) 
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For example, for neutron decay, the lifetime of neutrons T^, has the 
distribution function 

*(T d > = 1 - exp (-T d /r d ), ( 

where r, - 15,3 minutes is the neutron mean lifetime, 
a 


T d = -T d iri ( 1 _ V- ( A3 > 

Given the height h., the distribution function of h is the 

J J hi 

probability for the neutron to be scattered between h^ and h i.e. 

$(h j+l ) = 1 - es P ^-^ s (Ej) rN(hj) - N(h j+1 )]/u ] . (A4) 

Here CT g (E) is the cross section for elastic scattering at energy E, and 
N(h) is the columnar density of the solar atmosphere above the height h. 
From equations (Al) and ^A4) we get that 

N(h j+1 ) = N( V + ^j/ CT s ( E j)^ ^(l-R sJ ). (A5) 

By inverting N(h^_^) we obtain hj + ^ ; we note, however, that if 

and N(h^) <-^/a(Ej ) fn(l-R s ) no solution can be found. Physically 
this means that the neutron escapes after the scattering. 

Having found h^ +1 , the time t ^ is simply 

t j+i = c j + (h j+r h j )/(M j v j ) ( A6 > 

where v ^ is the velocity at energy E . The neutron will decay between 

t. and t if 
J J+l 

C j < T d < C j+l* ( A7 ) 
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In order to determine whether the neutron will be captured on 
protons between t^ and t^ +1 we proceed as follows. Let h cH be the 
height at which a neutron is captured after j elastic scatterings. 
Applying the same considerations as for equation (A4) , the distribution 


function of h^ is 


4(h cH ) - 1- 


j i exp r-o c (E k _ 1 )(N(h, t _ l )-N(h k ))/u t _ 1 ]|x 
X expr~a c (E j) (N(hj) - 


(A8) 


where cr = a (H) is the capture cross section on hydrogen as a function 
c c 


of neutron energy* Using equations (A5) and (A8) , we get that 
N(h cH ) = N(hj) + Mj/CT c (Ej) r ^ n < 1_R c H ) '^ =1 (E^) £n(1 " R s,k-l 


)], (A9) 


where R is another random numberjand the random numbers R g ^ are 
those which have been used to calculate h^. Equation (A9) can be 
inverted j the neutron will be captured on a proton between t^ and t^^^ 


if 


h - h 

< 

h - - h. 

cH j 


J+l J 


(A10) 


Simlarly, the neutron will be captured on a He between t and tj^ if 


h cHe 3 ' h j 


h j+l - h j 


(All) 


Here h tls obtained in the same manner as h , with the exception that 
cHe ^ cn 

3 3 3 

a (H) is replaced by (He /H) CT (He ) , where He /H is the ratio of the 
c v c 

He 3 to H in the photosphere. Once a neutron is captured by a proton. 
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we record the time and depth of capture « The former is used to 
calculate the time profile of the burst and the latter is needed to 
evaluate the effects of Compton scattering on the photons. In order to 
obtain the energy and direction of the neutron after elastic scattering 
at fit's t find the temperature at this height , and then use a 

Maxwell -Boltzman distribution to find the velocity components of the 
interacting protons (random numbers corresponding to gaussian distri- 
bution can be generated directly on a computer) * We then assume that 
the directional cosine of the neutron is uniformly distributed in the 
center of mass system so that we can choose a random number R 

np 

corresponding to this directional cosine. From the conservation of 
energy and momentum we then get the energy and direction of the scattered 
neutron in the frame of the Sun, 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Probabilities for neutron escape, decay, capture on protons 
3 

and loss on He in the solar atmosphere (solid lines) ; and 
photon yields per neutron (dashed lines) . The parameter 0 
is the angle between the earth-sun line and the heliocentric 
radial direction through the flare. The ratio He /H is the 
photospheric helium 3 abundance, and E^ is the energy of the 
neutrons . 

Probabilities for neutron escape, decay, and capture on 
protons in the solar atmosphere (solid lines) and photon 
yields per neutron (dashed lines) . 

Time profiles of the 2.2 MeV line intensity at Earth 

multiplied by 4 ttR , where R = 1 A.U., produced by instantaneously 

released monoenergetic neutrons at time t . The rest of the 

o 

parameters are the same as in Figures 1 and 2* 

Time profiles of the 2,2 MeV line intensity at Earth 
2 

multiplied by 4ttR 3 where R — 1 A,U., produced by instantaneously 

released monoenergetic neutrons at time t . The rest of the 

o 

parameters are the same as in Figures 1 and 2, 
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